SCIENTIFIC 

REPORTS 




OPEN 



SUBJECT AREAS: 

SOLAR ENERGY AND 
PHOTOVOLTAIC 
TECHNOLOGY 



Exchanging Ohmic Losses in 
Metamaterial Absorbers with Usefu 
Optical Absorption for Photovoltaics 



METAMATERIALS Ankit Vora', Jephias Gwamuri^, Nezih Pala^, Anand Kulkarni', Joshua M. Pearce''^ & Durdu O. Guney' 



Received 
28 October 2013 

Accepted 
10 April 2014 

Published 
9 May 2014 



Correspondence and 
requests for materials 
should be addressed to 
D.O.G. (dguney@mtu. 

edu) 



' Department of Electrical and Computer Engineering, Michigan Technological University, Houghton, Ml 4993 1 , USA, ^Department 
of Materials Science and Engineering, Michigan Technological University, Houghton, Ml 49931 , USA, ''Department of Electrical 
and Computer Engineering, Florida International University, Miami, FL 33 1 74 USA. 

Using metamaterial absorbers, we have shown that metallic layers in the absorbers do not necessarily 
constitute undesired resistive heating problem for photovoltaics. Tailoring the geometric skin depth of 
metals and employing the natural bulk absorbance characteristics of the semiconductors in those absorbers 
can enable the exchange of undesired resistive losses with the useful optical absorbance in the active 
semiconductors. Thus, Ohmic loss dominated metamaterial absorbers can be converted into photovoltaic 
near-perfect absorbers with the advantage of harvesting the fuU potential of light management offered by the 
metamaterial absorbers. Based on experimental permittivity data for indium gallium nitride, we have shown 
that between 75%-95% absorbance can be achieved in the semiconductor layers of the converted 
metamaterial absorbers. Besides other metamaterial and plasmonic devices, our results may also apply to 
photodectors and other metal or semiconductor based optical devices where resistive losses and power 
consumption are important pertaining to the device performance. 



Improvements in optical enhancement for solar photovoltaic (PV) devices using conventional natural materi- 
als has begun to show diminishing returns, yet metamaterials, which are rationally designed geometries of 
optical materials that can be tuned to respond to any region of the electromagnetic spectrum offer an 
opportunity to continue to improve solar device performance'. Historically, periodic structures with sub-wave- 
length features are used to construct a metamaterial for a given application. Designing and constructing a material 
in this way a metamaterial possesses optical properties that are not observed in their constituent materials 
enabling the properties to be determined from structure rather than merely composition. The advent of meta- 
materials has provided PV device designers among others with unprecedented flexibility in manipulating light 
and producing new functionalities^. Metamaterials have already been proposed for super lenses that allow sub- 
wavelength resolution beyond the diffraction limit, and electromagnetic cloaks, which promise real physical 
invisibUity-*. 

Another approach to optical enhancement of PV is to consider plasmonics, which is a rapidly growing field for 
the application of surface plasmons to device performance improvement. Surface-plasmon-based devices already 
proposed include not only PV but also light-emitting devices, data storage, biosensors, nano-imaging, wave- 
guides, perfect absorbers, and others^"". 

Surface plasmons are collective oscillations of surface electrons whereas surface plasmon polariton (SPP) 
describes a coupled state between a surface plasmon and a photon'. On the other hand, localized surface plasmons 
(LSP) are confined to bounded geometries such as metallic nanoparticles or nanostrips of various topologies'""'^. 
The fluctuations of surface charge results in highly localized and significantly enhanced electromagnetic fields in 
the vicinity of metallic surfaces. 

Typically, surface plasmon resonances exhibit a strong relationship to the size, shape and the dielectric 
properties of the surrounding medium. The resonances of noble metals are mostly in the visible or infrared 
region of the electromagnetic spectrum, which is the range of interest for PV applications'^. 

For the application, we are considering here photons arriving at the metal surface produce surface waves in the 
form of SPP along the metal-semiconductor interface of the top layer of a PV device. This occurs when the 
photons interact with the collective oscillations of free electrons in the metal of the absorber. The metallic 
nanostructures have the ability to maintain SPPs, which provides electromagnetic field confinement and 
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enhancement. This phenomenon has already attracted substantial 
attention in the PV scientific community, but has yet to be commer- 
cialized and scaled. 

Perhaps even more promising it is possible to construct plasmonic 
metamaterial "perfect absorbers" with such surface waves. 
Metamaterial absorbers are extremely flexible and have already been 
designed with broadband""", polarization-independent'^"^', and 
wide-angle'*'"'"'"'™'^^ optical absorption for a host of applications. 
Both theoretical work through device simulations'^''"'^^ and physical 
fabrication'"''"''"'""^''^' of ultrathin metamaterial absorbers with 
these features have been realized. The theoretical limit of nearly 
100% absorbance at the tunable resonant wavelength can be obtained 
with proposed designs of ultrathin, wide-angle perfect absorber 
structures for infrared light^^. Very recently, a study showed an aver- 
age measured absorption of 71% from 400-700 nm, which repre- 
sents a majority of the energy in the solar spectrum, with an ultrathin 
(260 nm), broadband and polarization independent plasmonic 
super absorber". In addition, the simulations of Aydin et al., indi- 
cated that the absorption levels could be increased to 85%". 

AU of these features are critical to maximize the efficiency of PV 
devices, yet are generally lacking in schemes to optically enhance 
solar cell design. Even small improvements in optical enhancement 
in specific cases can have a large effect. For example, improving 
(enlarging) the acceptance angle for high absorption for a PV device 
may eliminate the need for mechanical tracking of a PV system, 
which can have a substantial impact on both the first cost as well 
as the operating and maintenance costs. 

Complex perfect absorbers have more degrees of freedom in terms 
of impedance-match, polarization-independence, and wide angular 
reception, which make them superior systems for managing light 
than simple plasmonic nanoparticle approaches. Thus perfect 
absorbers are better candidates for enhancing solar PV conversion 
efficiency. 

Despite this promise, perfect absorbers have not been studied 
systematically in the literature for solar photovoltaic applications'. 
In addition, in the few cases where the potential was probed the 
device design limited the practicality of the approach. For example, 
Aydin, et al.'s "super absorber"", uses a lossless spacer where almost 
aU the power covering the solar spectrum is absorbed by the metal. 
This results in significant Ohmic loss via heating, which loses both 
the power of the photons but also decreases the efficiency of the next 
electrical conversion because of resultant higher cell operating tem- 
perature. Instead, the approach taken in this paper is to change the 
absorption localization from metal to the semiconductor layer by 
replacing the dielectric spacer with a "semiconductor spacer" in 
the perfect (or "super") absorbers, which can then be used to enhance 
thin-film PV device performance. 

The geometric skin depth (i.e., depth at which the current density 
is reduced to 1/e of its value at the surface of the conducting material) 
in a structured conductor such as metallic nanostructures (e.g., split- 
ring resonators, strips, etc.) is much larger than in the bulk form of 
the same conductor". However, in the limit the nanostructure 
approaches bulk conductor, the geometric skin depth is significantly 
reduced and approaches its bulk value. More importantly, this skin 
depth reduction can be exploited as an important resource to min- 
imize Ohmic losses in metallic nanostructures despite higher metal- 
lic filling ratio. We will call this technique the "Bulk Skin Depth 
Technique (BSDT)." Previously, we applied this technique to reduce 
Ohmic losses in metamaterials^* and to design three-dimensionally 
isotropic negative index metamateriaP' (see Fig. SI, Supplementary 
information). Now, in this paper we apply the BSDT to perfect 
absorbers to turn them into functional solar cell devices and use 
the indium gallium nitride (In-^Gai.xN) PV device as an example. 

If we reduce the geometric skin depth in the metallic portions of 
the metamaterial absorbers (either perfect or near-perfect), we can 
push the resonant currents induced by solar radiation closer to the 



semiconductor layers in PV devices. This brings the following 
important advantages: (1) Optical absorption in the metallic layers 
(i.e., Ohmic loss) is shifted into semiconductors, which increases 
useful optical absorption for solar energy conversion. (2) Joule heat- 
ing is minimized. (3) Because the current is tightly confined to metal- 
lic surfaces facing the semiconductor layers, additional interconnects 
or contacts can be made possible^' without destroying or short-cir- 
cuiting the absorber (see Fig. SI, Supplementary information). 

In order for a plasmonic PV cell to reach the full potential of the 
technology the designed cell must maximize the absorption within 
the semiconductor region. Simultaneously, absorption must be mini- 
mized everywhere else including the metallic regions over the cell. 
This balancing act of absorption must be accomplished over the 
entire solar spectrum (normally AMI. 5, which is ASTM G-173) 
and ideally over all incidence angles. Previous work has shown that 
a BSDT" can be employed to minimize Ohmic losses in metamater- 
ial absorbers by guiding optical absorption to the semiconductor 
layers and away from the metal. In order for this technique to func- 
tion, it first requires identification of resonant modes of the meta- 
material absorber that contributes to optical absorption. Then the 
nanostructure of the metal must be geometrically customized as 
determined by the underlying resonant currents. 

The BSDT approach can be applied to different geometries (e.g. a 
square-grid structure made up of cross structures formed by perpen- 
dicularly placed metallic strips) or any other polarization-independ- 
ent metamaterial absorber design found in the literature"'"'^"''^'''^'. In 
addition this approach will work for both "super absorber"" as well 
as metamaterial absorbers with broadband response employing the 
concept of nanoantenna'*'". 

Results 

Applying the BSDT to a metamaterial perfect absorber and 
exchanging Ohmic losses. Below we show how to apply the BSDT 
to a simple perfect absorber to illustrate the concept. Application to 
other perfect absorbers is similar. In Fig. SI (see Supplementary 
information) we provide a different (earlier) example where the 
BSDT is applied to split-ring-resonators with resonant loop 
currents". Figs. l(a)-(c) show a two-dimensional cross section of 
the simple metamaterial absorber, induced resonant currents, and 
power loss density, respectively. The absorber consists of gold (red) 
layers and a lossless spacer (blue). 

In most metallic nanostructures operating around optical frequen- 
cies the geometric skin depth is larger than both the bulk skin depth 
and the dimensions of the cross-sectional area of the metal through 
which induced currents flow. This usually results in larger Ohmic 
losses at optical frequencies compared to lower frequencies due to 
deeper penetration of the photons relative to structure dimensions". 
Therefore, it is desired to reduce the geometric skin depth at optical 
frequencies to reduce the Ohmic losses. The BSDT offers an effective 
strategy for the geometric tailoring of the skin depth. In order to 
apply the BSDT we first need to identify the underlying resonant 
mode of the nanostructure and then modify the geometry such that, 
for induced resonant currents, the metallic surfaces of the nanos- 
tructure behave effectively similar to the surfaces of semi-infinite 
bulk metals with bulk skin depth. Below, we explain how the 
BSDT is applied to the metamaterial absorber sketched in Fig. 1(a). 
Fig. 1(b) shows the current density distribution for the identified 
plasmonic resonant magnetic dipole mode responsible for perfect 
absorbance. The dipole mode has antisymmetric current flow at 
the top metal strip and bottom ground plate indicated by black 
arrows. The current density has a maximum at the metal-spacer 
interface, where the surface plasmons reside, and decays exponen- 
tially inside the metal with the increasing distance from the interface. 
Fig. 1 (f) shows that when fj = 20 nm the geometric skin depth for the 
strip is 38 nm which is larger than both bulk skin depth (i.e., 23 nm) 
and the cross-section of the strip through which the current flows. To 
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Figure 1 | (a) Two-dimensional cross-section of a unit cell of a simple absorber consisting of a lossless spacer (shown blue) sandwiched between a 
strip and a ground plate both made from gold (shown red) . The size of the unit cell along the direction of periodicity is p = 350 nm. The width of the strip 
is w = 250 nm. The strip and the ground plate have the same thickness of tj = = 20 nm. The spacer thickness is d= 14 nm. The electric permittivity of 
the spacer is E = 2.25. (b) Resultant magnetic dipole oscillating at 150 THz. The antiparallel current density (A/m^) indicated by arrows give rise to optical 
absorption in metallic regions as evidenced by (c) the power loss density (W/m'). (d) Absorbance at different regions of the perfect absorber with the same 
design parameters as (a)-(c) except tg = 60 nm to prevent transmittance. Total absorbance is almost 100% shared between the strip (—70%) and the 
ground plate (—30%), since lossless spacer is assumed, (e) Using BSDT total absorbance in the metal is reduced to —70% shared equally between the strip 
and the ground plate. BSDT is applied by increasing metallic thicknesses to 1^= tg = 85 nm. (f) Geometric skin depths for the strip and the ground plate 
for different strip thicknesses. The ground plate has a fixed thickness of = 75 nm except the last data point where = 85 nm. Geometric skin depth 
converges to bulk skin depth with increasing metallic thicknesses, (g) Introducing a non-zero imaginary part to the permittivity of the original lossless 
spacer in (d)-(e) (e.g., £ = 2.25 + 0.2 j) leads to the change of absorption localization from metals to the spacer. This lossy spacer can be replaced with a 
semiconductor material for solar cell applications to exchange metallic Ohmic losses with useful optical absorption for carrier generation in the 
semiconductor regions. 



reduce the skin depth we need to tailor the geometry such that the 
finite metal surfaces, where the resonant currents are bounded to, 
approach to bulk metal surfaces so that the geometric skin depth 
approaches bulk skin depth (or even below). For the given structure 
in Fig. 1(a), without an important shift in the resonance frequency, 
this can only be achieved by increasing and tg. Although, at first, 
increasing the w might also be seen accommodating for the "bulk 
limit," we cannot change w since it is one of the critical parameters to 
tune the resonance frequency of the magnetic dipole as we will 
explain later when we discuss the scalability of the perfect absorber 
toward the red edge of the visible spectrum. 

Fig. 1(d) shows near 100% total absorbance achieved by the simple 
perfect absorber using the same design parameters as Figs. l(a)-(c) 
except the ground plate thickness = 60 nm to prevent transmit- 
tance. When the spacer is assumed to be a lossless dielectric, the total 



absorbed optical power is shared only between the top and bottom 
gold layers. Using BSDT, we reduce the total absorbance in the metals 
to about 70% as shown in Fig. 1(e). Here, BSDT is applied by making 
metallic layers thicker along the direction perpendicular to the flow 
of conduction currents; because this allows the metallic surfaces 
behave like semi-infinite bulk metal surfaces and reduce the geomet- 
ric skin depth for the induced resonant currents, hence reduce the 
Ohmic losses similar to Fig. SI (see Supplementary information). 

Orange line (diamond symbols) in Fig. 1(f) shows the geometric 
skin depth for the strip approaches to bulk skin depth as the thickness 
of the strip increases. At the same time, the peak absorbance in the 
strip drops from about 70% at tg = 20 nm [see Fig. 1(d)] to about 35% 
at = 85 nm [see Fig. 1(e)]. The ground plate thickness is fixed at tg 
= 75 nm during the process except at the last data point where tg = 
85 nm. This does not affect the skin depth trends for the strip or the 
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ground plate as they have already converged to the bulk skin depth. 
The blue horizontal line (square symbols) corresponds to the geo- 
metric skin depth for the ground plate which is almost constant and 
about the same as bulk skin depth due to sufficiently large thickness. 
We should note that our observation in Fig. 1(f) is consistent with 
Fig. 7 in Ref. 24. 

The most of the unabsorbed power in Fig. 1(e) is reflected due to 
impedance mismatch and the localization of some absorbed power is 
internally changed from the strip to the ground plate since BSDT is 
mainly applied to the strip here. Note that the strip thickness 
increases 65 nm during the BSDT process while the ground plate 
thickness increases only 10 nm. In general, if we apply the BSDT to 
only a subset of metallic components, internal change of absorption 
localization from those metallic components to other metallic com- 
ponents occur as we observe here. This internal redistribution of 
power dissipation occurs at a lesser degree if BSDT leads to imped- 
ance mismatch. 

We can take advantage of this impedance mismatch in the next 
stage by replacing the lossless spacer with a lossy spacer (such as a 
semiconductor). For example, if we introduce a non-zero imaginary 
part to the permittivity of the original lossless spacer [see Fig. 1(g)], 
the reflected power can be acquired back into the lossy spacer, hence 
not only the Ohmic losses are reduced but also the optical absorption 
is restored back to nearly 100%, maximized in the lossy spacer [com- 
pare with Figs. l(d)-(e)]. The lossy spacer can improve the imped- 
ance mismatch which was created by BSDT. If the imaginary part of 
the electric permittivity of the lossy spacer is sufficiently large, almost 
aU the power originally absorbed by the metal can be shifted to the 
spacer as discussed in our paper. 

Metamaterial absorbers with IiixGai.xN spacers. To explore the 
applicability of the technique to the semiconductors used in PV 
technology, we chose InxGai.xN as the spacer because 1) it has a 
potential band gap engineering due to variable band gap from 0.7 
to 3.4 eV created by changing the indium concentration that covers 
nearly the entire solar spectrum, 2) it has a high absorption 
coefficient that would make it possible to utilize thin absorber 
layers most applicable to this method of optical enhancement, 3) 
recent microstructural engineering has been demonstrated with 
nanocolumns that offer optical enhancement, and 4) this provides 
the potential for several known photovoltaic device configurations 
and multi- junctions with theoretical efficiencies over 50%^""'". 

The experimental complex electric permittivity (e) data for 
Ino.54Gao.46N'' is given in Fig. S2 (see Supplementary information). 
We applied the BSDT at five specially selected frequencies. These 
frequencies along with the corresponding e for Ino.54Gao.46N are 
listed in Table 1. If it is desired, it should be possible to apply the 
technique to any other intermediate frequencies between these five 
representative frequencies. Because these frequencies are selected 
carefully by considering the minima, maxima and profile of the 
complex permittivity data as described below, the absorbance results 



for the intermediate frequencies should be interpolations of the 
results obtained from the selected frequencies. 

The complex permittivity data in Fig. S2 (see Supplementary 
information) starts from 145 THz, at which the permittivity is ei 
= 0.86 + 1.23 j. This is selected as the first frequency point. 
Above this point, the 9(e) monotonicaUy decreases and reaches 
global minimum at 387 THz while the Jf (e) manifests a monotonic 
increase. Thus, the permittivity at this second selected frequency 
point becomes 82 = 4.59 -I- 0.45 j. The third selected frequency point 
is reached at 695 THz where the permittivity is 83 = 5.67 -I- 1.61 j. 
This is the frequency point where 5R(£) reaches a global maximum 
and the 9(e) still tends to increase. At the fourth selected frequency 
point, 1000 THz, the 3ff(e) reaches a minimum while the 9(e) 
reaches global maximum. The permittivity at this point becomes 84 
= 5.06 + 2.74 j. The experimental data ends at 1090 THz where the 
permittivity becomes 85 = 5.47 + 2.60 j. This is selected as the fifth 
special frequency point. 

Table 1 also lists geometric parameters for five metamaterial 
absorbers tuned to operate around the five respective selected fre- 
quencies (referred as Points 1-5) listed in the table. The geometry for 
the assumed metamaterial absorbers and the definitions of the geo- 
metric parameters are described in Figs. l(a)-(c). The spacers used in 
the absorbers are initially chosen to be lossless and their permittivity 
is equal to the 5R(e) at the respective selected frequency. Below we 
show that these metamaterial absorbers can be turned into near- 
perfect "solar absorbers" where most of the metallic absorption is 
shifted to semiconductor layers (i.e., InxGai.xN) by using BSDT fol- 
lowed by replacing the lossless spacers with the semiconductor 
spacers. 

Exchanging Ohmic losses at the infrared spectrum. First, we apply 
the BSDT at around Point 1 to shift the optical absorption in the 
metallic regions to the semiconductor layer. We start from a perfect 
absorber which has a peak absorbance at about 180 THz [see 
Fig. 2(a)], relatively close to Point 1. The thickness of the ground 
plate was chosen as 60 nm to prevent transmittance. The other 
geometric parameters for this perfect absorber are given in Table 1. 
Figs. 2(b) and (c) show the magnetic field and current density 
distribution, respectively, for the resonant magnetic dipole mode 
leading to the perfect absorption. Relatively simple field profile of 
the resonant mode renders the application of BSDT fairly 
straightforward. Because the surface plasmons are mainly confined 
to metal-spacer interfaces, we can push the surface plasmons further 
to the interface by increasing the thicknesses of both the strip and the 
ground plate. This wiU reduce the skin depth and prevent further 
penetration of the electromagnetic fields into the metallic regions. 
When the strip thickness reaches from its original value of 15 nm to 
the same thickness as the ground plate, the photons cannot penetrate 
deep into the strip anymore and hence the peak absorbance in the 
strip reduces from about 75% to 30% [see Fig. 2(d)]. The photons 
which are not absorbed anymore by the strip are either reflected or 
added to the absorbance in the ground plate. The overall peak 
absorbance of about 60% is shared almost equally by the strip and 
the ground plate at the same thicknesses of 60 nm. The overall final 
peak absorbance can be reduced further down to about 48%, again 
shared almost equally by the strip and the ground plate when they 
both reach 245 nm thickness [see Fig. 2(e)]. The unabsorbed 
photons in the process are reflected. 

Having minimized the metallic absorbance in the perfect absorber 
by BSDT, in the next stage we gradually restore the 9(e) of 
Ino54Gao.46N in Fig. 3. This leads to the total absorbance in the 
metallic regions being further minimized to below 15% at the first 
targeted frequency 145 THz. The photons expelled from the metallic 
regions as well as photons reflected originally due to impedance 
mismatch with the increasing metal thicknesses during the applica- 
tion of BSDT are now shifted to and absorbed in the semiconductor 



Table 1 | Real and imaginary parts of the complex permittivity for 
the lno.54Gao.46N data in Fig. S2 (see Supplementary information) 
at five selected frequency points and geometric parameters for the 
five perfect (or near-perfect) absorbers operating near those five 
respective selected frequencies 



Data Point 


Frequency 
(THz) 


R(£) 


3(E) 


p(nm) 


w [nm] 


h (nm) 


tg [nm] 


d (nm) 


Point 1 


145 


0.86 


1.23 


400 


315 


15 
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14 


Point 2 
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4.59 


0.45 


150 


55 
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Point 3 
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5.67 


1.61 
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2.74 
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Figure 2 | (a) Absorbance in different parts of the perfect absorber operating in the neighborhood of Point 1 in Table 1. (b) Magnetic field (A/m) 
component perpendicular to the cross section (i.e., defined as H^) and (c) distribution of the current density (A/m^) parallel to the metal-spacer interface 
(i.e., defined as /J at the peak absorbance. Red and blue colors correspond to the current flowing to the right and left, respectively, (d) Absorbance in 
difl^erent parts of the absorber when the strip thickness increases to = 60 nm and (e) when both the strip and the ground plate have 1,= tg= 245 nm 
thickness. All the other parameters are kept fixed and the same as in Table 1. 



layer (see Fig. 3). The total peak absorbance in the semiconductor 
layer at the target frequency exceeds 85% and the overall perfect 
absorption capability of the absorber is restored when the 9(e) takes 
its true experimental value of 1.23. In contrast, the maximum total 
absorbance in plain Ino.54Gao.46N layer for the same thickness (i.e., d 
= 14 nm) is only about 6% based on available experimental complex 
permittivity data in Fig. S2 (see Supplementary Information). Thus, 
we convert the Ohmic loss dominated perfect absorber to "photovol- 
taic near-perfect absorber." In this stage, the resonant surface plas- 
mon mode displayed in Figs. 2(b) and (c) again plays an important 
role on changing absorption localization from the metallic regions to 
the semiconductor layer. Since the surface plasmons reside at the 
metal- semiconductor interface, the accompanying photons can be 
easily absorbed near the interface by the semiconductor layer. This 
process of shifting absorption localization to the semiconductor 
layers can be achieved over the entire infrared spectrum under the 
same strategy above using the same magnetic dipole mode. 

Below, we study the role of BSDT in more detail to better under- 
stand the transition from perfect metallic absorbance in Fig. 2(a) to 
near-perfect semiconductor absorbance in Fig. 3(d). Fig. 4(a) shows 
absorbance in different parts of the absorber and normalized skin 
depth (i.e., normalized with respect to its maximum value of 53 nm) 
versus in the same graph. The skin depth is maximized at small 
values and decays exponentially with increasing and finally con- 
verges to bulk skin depth values when the becomes sufficiently 
large. The change in the skin depth has important impact on the 
surface impedance'^'^^ of the metamaterial absorber. Below, we show 
that tailoring the geometric skin depth can be exploited as an import- 



ant and convenient resource to control the surface impedance and 
absorbance in different parts of the absorbers. 

As shown in Fig. 2(a) the structure behaves as an impedance 
matched perfect absorber with the lossless spacer when fj = 15 nm. 
Because, when = 15 nm, the skin depth is 53 nm (i.e., larger than 
the tj). This allows photons to easily penetrate through metal and be 
completely absorbed as shown in Fig. 2(a). However, when the 
becomes larger, the metal strip approximates bulk metal and thus 
the skin depth for the resonant magnetic dipole mode decreases and 
approaches to bulk skin depth as shown in Fig. 4(a). Since the skin 
depth is reduced, photons cannot penetrate deeper into the metal. For 
example, when = 35 nm, the skin depth is reduced to 32 nm (i.e., 
less than the Q. This not only reduces the absorbance but also 
increases the impedance mismatch. Increasing the further above 
60 nm results in converging skin depth. Since the skin depth is con- 
verged, both the absorbance and impedance mismatch also tend to 
converge as can be seen from absorbance plots in Fig. 2. In the fol- 
lowing we explain that the impedance mismatch achieved by tailoring 
the skin depth can be used to control where the absorbance is localized. 

The absorbance curves in Fig. 4(a) are obtained by replacing the 
lossless spacer in Fig. 2 with the lossy spacer such that 
Ji(e) = 5R(ei ) = 0.86 and 3(e) is assumed to have four different values 
which are 0.1, 0.3, 0.6, and Q?(ei) = 1.23, to illustrate the interplay 
between the skin depth, imaginary part of the permittivity of the lossy 
spacer, and impedance matching. Att^ = 15 nm we observe that the 
peak absorbance in the lossy spacer is decreased with increasing 9(e), 
because the absorber [see Fig. 2(a)] at = 15 nm with lossless spacer 
was an impedance matched perfect absorber. When we increase the 
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Figure 3 | Absorbance in different parts of the absorber achieved by the BSDT process described in Fig. 2 followed by the gradual incorporation of 
the true complex permittivity such that (a) Q(£) =0.1, (b) Q(£) =0.4, (c) Q{s) =0.8, and (d) Q(£) = 9(£i) = 1.23. The total absorbance in plain 
Ino.54Gao layer of the same thickness (i.e., d = 14 nm) is also shown for comparison using the experimental complex permittivity data in Fig. S2 (see 
Supplementary information) (e) i-f, (A/m), (f) }^ (A/m^), and (g) the power loss density (W/m') distribution for the near-perfect absorber in (d) for the 
first selected frequency of 145 THz in the Ino.54Gao.46N data (i.e.. Point 1). Significant portion of the power (i.e., over 85%) is absorbed inside the 
Ino.54Gao.46N layer. 



imaginary part of the lossy spacer, the structure becomes increasingly 
impedance mismatched. Therefore, the absorbance in both lossy 
spacer and metals decreases with increasing Fig. 4(b) shows 

surface impedance for Q'(£i ) = 1 .23 and normalized skin depth versus 
ts in the same graph. Notice that in the region where the skin depth is 
high, the structure is impedance mismatched and starts to become 
impedance matched as the skin depth converges. (See Fig. S3, 
Supplementary information, to analyze the impedance matching for 
different Q?(£) values based on total absorbance plots.) 

In Fig. 4(a), for = 155 nm, we can rank the absorbance in the 

lossy spacer for different imaginary parts as Aq \ > Aq \ >Aq \>A [^\^ 



where the superscript and subscript stands for the spacer and 3(e), 
respectively. On the other hand, at = 35 nm, we observe that the 
absorbance starts to drop after a certain value of the 9'(£), that is 

^03>^o'i but A[%<A''gl<A''ol. This means that the impedance 
mismatch produced by the reduced skin depth at tj = 35 nm is only 
to the extent that 9(e) = 0.3 and other larger imaginary values over- 
compensate the impedance mismatch and thus result in smaller 
absorbance. Therefore, the optimum value of 9(e) should be some- 
where between 0.3 and 0.6 to maximize the absorbance in the lossy 
spacer. This optimum value of the 9(e) is the one which fully com- 
pensates the impedance mismatch opened by reduced skin depth at 
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Strip Thickness (nm) 

Figure 4 | (a) Absorbance and normalized skin depth (NSD) versus tj. The purple curves show the peak absorbance in the spacer region. The 
green curves show the total absorbance in the metallic regions. Orange curve shows the skin depth normalized with respect to its maximum value of 
53 nm at = 15 nm. aJ'' stands for absorbance in spacer (i = s) or metaUic regions (j = m) for different Q(£) values j of the spacer. NSD is calculated at 
Q(£) = 0. 5?(£) = 3?(£i) =0.86 and = 245 nm. (b) Effective surface impedance of the metamaterial absorber and NSD [repeated from (a)] versus tj 
calculated near the total absorbance peak for £i = 0.86 + 1.23 j and = 245 nm. (c) The percentage of total absorbance in different regions of the 
metamaterial absorber corresponding to (a) and NSD versus t^. stands for the percentage of total absorbance in region i for different Q(£) values ; of the 
spacer. 



= 35 nm. Similary, if we consider = 150 nm, the resultant total 
impedance mismatch should be compensated by an Q?(£) between 0.6 
and ^(ei) = 1.23. Since the impedance mismatch increases with t^, 
the optimum value of Sy(£) also increases. We noticed that our obser- 
vations here are consistent with a recent study by Hu, et al.^' where 
they study the role of Q?(£) of the spacer on the impedance match and 
absorption. What we add here is the effect of changing skin depth 
which is not considered in Ref. 33. 

Another important observation in Fig. 4(a) is that the curves for 
9 (a) corresponding to the absorbance in the lossy spacer first 
increase and then decrease except for 9(e) = 9'(ei) = 1.23. The 
increase is rapid and occurs in the regime where the skin depth 
and impedance change rapidly whereas the decrease is more gradual 
[no decrease at all for ^(fi) = 9(£i) = 1.23] and occurs in the regime 
where the skin depth and impedance tend to converge. The peaks for 
different imaginary parts show that the specific value of 9'(g) cannot 
compensate the impedance mismatch if is increased further 
beyond the peak, because this results in smaller skin depth and/or 



higher impedance mismatch which requires larger 3(e) for full com- 
pensation of the resultant impedance mismatch. The peak for 
9(e) = Sy(ei) = 1.23 occurs at = 245 nm and corresponds to the 
peak in Fig. 3(d). 

Regarding metallic losses in Fig. 4(a) we observe in general that the 
metallic losses decrease with increasing 9(e) . However, the change in 
the skin depth and impedance has also an impact on the metallic 
losses. For example, consider 9(e) = 0.1. At = 15 nm the structure 
with the lossless spacer was an impedance matched perfect absorber 
and had the maximum skin depth. Under this condition. Fig. 4(a) 
shows that the incorporation of any lossy spacer with the given 
relatively large 9(e) values increases the impedance mismatch. The 
smaller the 9(e), the smaller is the impedance mismatch. Among the 
given 9(e) values, 9(e) = 0. 1 is the smallest. Therefore, the reduction 
in metallic losses is the smallest for 9(e) = 0.1, in turn followed by 
9(e) =0.3, 9(e) =0.6, and 9(e) =9(ei) = 1.23. For 9(e) =0.1, if we 
increase the t^, the metallic losses gradually first decrease and then 
increase. The initial decrease happens, because in this regime the 
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metallic losses rapidly decrease as a result of rapid decrease in the 
skin depth [see Figs. 2(a) and (b)]. Additionally, increase in the 
impedance mismatch with increasing in this regime also plays a 
role in reducing the metallic losses. However, at larger values the 
contribution of the 9(g) = 0. 1 to compensate the resultant mismatch 
decreases. Therefore, the metallic losses at these values cannot be 
reduced further below their original values when the structure was 
having lossless spacer. This also applies to larger 9(e) values, which 
also exhibit increasing metallic losses with increasing t^. On the other 
hand, in contrast with 9(e) =0.1, the larger 9(e) values show rela- 
tively flat and low metallic absorbance at small values. Because, 
initially at = 15 nm, the structure has substantially large imped- 
ance mismatch for all the 9(e) values. Therefore, the absorbance in 
both lossy spacer and metals starts from the lowest values at = 
15 nm. Increasing the reduces the metallic losses due to reducing 
skin depth. However, simultaneously improved impedance match 
increases the metallic losses. This is why the skin depth balances 
the impedance match in the flat region for metallic losses in Fig. 4(a). 

Finally, Fig. 4(c) shows how the percentage of total absorbance 
changes in the lossy spacer and metals with t^. Att^ = 15 nm, if the 
structure would have lossless spacer, the metals would have 100% of 
the total absorbance whereas the spacer would have 0% [i.e., perfect 
absorber. Fig. 2(a)]. When we replace the lossless spacer with the 
lossy spacer with 9(e) =0.1, we observe that the percentage of total 
absorbance in the lossy spacer (i.e., 56%) is higher than the metal (i.e., 
44%). Higher percentages are obtained in the lossy spacer for larger 
9(e) despite smaller absolute absorbance due to decreasing imped- 
ance match at = 15 nm. As a result of overall interaction between 
lossy spacer and metals summarized in Fig. 4(a), we observe an 
interesting exchange of absorbance localization between the lossy 
spacer and metals as we increase the in Fig. 4(c). The percentage 
of total absorbance in the lossy spacer (metals) first increases 
(decreases) in the region where the skin depth rapidly changes and 
then decreases (increases) in the region where the skin depth 
converges. 

Scaling the Ohmic loss exchange toward the red edge of the visible 
spectrum. To demonstrate the scalability of the photovoltaic 
absorber to Point 2 in Table 1, we start from a perfect absorber 
which is now tuned to a frequency around Point 2 by modifying 
the geometric parameters of the perfect absorber in Fig. 2(a). 
Considering the underlying resonant magnetic dipole mode 
displayed in Figs. 2(b) and (c), we can crudely estimate the 

resonance frequency from j v/lC where L and C are 

equivalent circuit inductance and capacitance, respectively (More 
accurate expression requires detailed consideration of the mode 
distribution, but this does not change our conclusions drawn 
here). Because the 5R(e2)=4.59 is substantially large compared to 
5R(ei) = 0.86, scaling the resonant magnetic dipole mode to the 
neighborhood of Point 2 requires a sufficiently narrow strip width 
of 55 nm to compensate increasing capacitance. The other geometric 
parameters mainly serve for optimization of the perfect absorber. 
The absorbance spectrum of the resultant perfect absorber is 
shown in Fig. 5(a). The trend we observe with the increasing metal 
thicknesses is similar to around Point 1 [see Figs. 5(a)-(d)]. This is 
due to the same electromagnetic mode tailored using BSDT. When 
the strip thickness gradually changes from 10 nm to 135 nm, the 
peak absorbance for the strip drops from above 75% to below 40%, 
while the absorbance in the ground plate, first, like before, increases 
to about 30% to accommodate photons not absorbed by the strip and 
then drops back to a level slightly above 20% with increasing 
thickness. Thus BSDT gives about 40% overall decrease in the 
metallic absorbance with lossless spacer. The magnetic field and 
current density distribution for the final absorber at the peak 
absorbance is shown in Fig. S4 (see Supplementary information). 



Despite thick metals the fields are tightly confined to the metal- 
semiconductor interfaces. 

Restoring gradually the true value of 9(e2) =0.45 shifts most of 
the remaining metallic loss to the absorption in the semiconductor 
layer [see Figs. 5(e)-(g)]. The overall peak absorbance in the semi- 
conductor layer reaches above 70% while the total metallic absor- 
bance drop to below 30%. The maximum total absorbance in plain 
Ino.54Gao.46N layer with the same thickness (i.e., d = 12 nm) is only 
about 5% for the given frequency spectrum. Power loss density dis- 
tribution in Fig. S5 (see Supplementary information) shows signifi- 
cant concentration of optical power in the semiconductor layer. 

In comparison with Point 1, photovoltaic absorber at Point 2 has 
about 15% less peak absorbance in the semiconductor. This is mainly 
due to relatively small 9(e2) =0.45 compared to 9(ei) = 1.23. Peak 
absorbance in the semiconductor layer should vary from about 85% 
to 70% between Point 1 and Point 2, respectively, due to the smooth 
variation of the complex permittivity between Point 1 and Point 2. 

Below, similar to Point 1, we study the role of BSDT on the trans- 
ition from perfect metallic absorbance in Fig. 5(a) to near-perfect 
semiconductor absorbance in Fig. 5(g). In Fig. 6(a) we study the 
absorbance characteristics of the structure around Point 2. 
Fig. 6(a) is obtained by replacing the lossless spacer with the lossy 
spacer similar to Fig. 4(a) such that 5R(e) = ^{£2} =4.59 and 9(e) is 
assumed to have two different values which are 9(e) =0.1 and 
9(e) = 9(e2) =0.45. At = 10 nm, the structure with the lossless 
spacer is an impedance matched perfect absorber [see Fig. 5(a)]. 
Here, if the lossless spacer is replaced with a lossy spacer with suffi- 
ciently small 9(e) (e.g., 0.1) the lossy spacer cannot open an imped- 
ance mismatch. The total absorbance is shared between the metals 
and lossy spacer. In this case, when 9(e) =0.1 and = 10 nm, the 
metallic absorbance drops to 69% from 100% and the absorbance in 
the spacer increases to 28% from 0%. However, if the 9(e) is suffi- 
ciently large [e.g., 9(e2) =0.45], the impedance mismatch occurs, 
metallic losses are further reduced to 27% and absorbance in the 
lossy spacer increases to 48%. Thus the total absorbance drops to 
75% due to the impedance mismatch. If we would further increase 
the9(e), the situation would be similar to fj = 15 nmin Fig. 4(a) [i.e., 
the larger is the 9(e), the higher is the impedance mismatch, the 
lower is the absorbance in both metals and lossy spacer.)] When we 
increase the t^, the impedance mismatch for the structure with loss- 
less spacer increases. Larger 9(e) compensates the opened imped- 
ance mismatch better than small 9(e) through absorbance in the 
lossy spacer. This is similar to our observation in Fig. 4(a). 
Fig. 6(a) also provides an interesting observation which does not 
occur in Fig. 4(a). At sufficiently small 9(e) the impedance com- 
pensation is mainly achieved by absorbance in metals, because the 
9(e) is not sufficiently large to reduce the metallic losses. However, at 
larger 9(e) values, the impedance compensation is achieved by 
absorbance in the lossy spacer. Another difference between 
Fig. 4(a) and Fig. 6(a) is that the absorbance curves tend to converge 
with the converging skin depth. Fig. S6 (see Supplementary informa- 
tion) shows that the trends for the percentage of total absorbance for 
both metals and lossy spacer in the regime where the skin depth 
rapidly changes are similar to Fig. 4(c), but the percentages converge 
with the converging skin depth consistent with the absorbance 
curves at this frequency range around Point 2. 

Briefly, Figs. 4(a) and 6(a) show that there are different regimes 
depending on the value of the 9(e) for a given magnetic dipole mode 
based absorber. If the 9(e) is sufficiently small, the incorporation of 
the lossy spacer does not change the impedance, but the metallic 
losses are not reduced much either [see, for example, = 10 nm 
in Fig. 6(a)]. If the 9(e) is gradually increased to the extent that the 
impedance mismatch is compensated, the metallic losses decrease 
and the absorbance in the lossy spacer increases [see, for example, 
= 145 nm in Fig. 4(a)] . If the 9(e) is increased beyond the point that 
the impedance mismatch is overcompensated, both metallic losses 
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Figure 5 | Absorbance in different parts of the perfect absorber (a) operating in the neighborhood of Point 2, (b) when the strip thickness in the 
absorber increases to = 60 nm, (c) when both the strip and the ground plate have = = 60 nm thickness, and (d) tj = = 135 nm thickness. All the 
other parameters are kept fixed and the same as in Table 1 . Absorbance in different parts of the absorber with the parameters same as in (d) except for (e) 
Si(E) =0.1, (f) Q(£) =0.2, and (g) Q(£) =Q(£2) = 0.45. The total absorbance in plain Ino54Gao.46N layer of the same thickness (i.e., d = 12 nm) is also 
shown for comparison. 



and absorbance in the lossy spacer is reduced [see again = 145 nm 
in Fig. 4(a)]. The ideal condition to maximize the absorbance in the 
lossy spacer occurs in the intermediate regime. The BSDT provides a 
simple guiding mechanism to achieve this ideal condition by con- 
trolling the impedance prudently in a desired fashion. 

Finally, Fig. 6(b) shows that the shift in the resonance peaks due to 
the change in geometric or material parameters does not have 
important effect on the observed distributions of the power dissipa- 
tion. To demonstrate this we kept the resonance frequency at 
378 THz [i.e., the resonance frequency when = 45 nm in 
Fig. 6(a)] for all the ts values by changing the original width (i.e., 
w = 55 nm) of the top metal strip. The maximum shifts in the 



resonance frequency and w from their original values occurred at 
= 15 nm. The shifts were 48 THz and 10.3 nm for the resonance 
frequency and tv, respectively. 

Scaling the Ohmic loss exchange to the violet edge of the visible 
spectrum using higher order plasmonic modes. Next, we scale the 
Ohmic loss exchange to Point 3, 695 THz (i.e., violet wavelength). 
Despite relatively large ^{£3) = 5.67, at this frequency, it is difficult to 
use the same magnetic dipole mode which requires vanishing strip 
width. Rather we can employ higher order plasmonic modes. 
Magnetic field distribution for one such mode is displayed in 
Fig. 7(a). This mode can be used to design a near-perfect absorber 
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Figure 6 | (a) Absorbance and NSD versus t^. The purple curves show the peak absorbance in the spacer region. The green curves show the total 
absorbance in the metallic regions. Orange curve shows the skin depth normalized with respect to its maximum value of 86 nm at = 10 nm. Aj'' stands 
for absorbance in region z for different SJ( a) values; of the spacer. NSD is calculated at Q(e) =0. ?R{e) =5R(£2) =4.59 and tg= 135 nm. (b) The absorbance 
[a!'' ] and the percentage of total absorbance [of] in different parts of the metamaterial absorber in (a) for Q(£) = Q(£2) = 0.45 versus t^ compared with 
the absorbance [Aj''*^' ] and the percentage of total absorbance [oj'"'^' ] of the metamaterial absorber with the constant resonance frequency at 378 THz. The 
curves for Q(£) = 0.1 were not shown for clarity. 



[i.e., above 95% peak absorbance, see Figs. 7(b)-(d)] in the 
neighborhood of Point 3 with a lossless spacer by tailoring the 
geometric parameters according to Table 1. Importantly, surface 
plasmons reside at the metal- semiconductor interfaces. This allows 
us to apply BSDT as before by simply increasing the metal 
thicknesses. Figs. 7(b)-(d) show total metallic peak absorbance 
drops down to about 40% by increasing the strip thickness from 
20 nm to 60 nm (i.e., reaching the same thickness as the ground 
plate). Here, in contrast with Points 1 and 2, first, the peak 
absorbance in the ground plate of fixed thickness decreases with 
the increasing strip thickness due to the increasing impedance 
mismatch at the corresponding frequencies. Second, the metals do 
not have to be very thick. When we incorporate the true value 
Q{si) — l.61 following the BSDT, the absorbance in the 
semiconductor layer at 695 THz exceeds 85% [see Figs. 7(e)-(g)]. 
In contrast, the total absorbance in plain lug 54Gao layer for the 
same thickness (i.e., d = 31 nm) is between the ranges of 15%-27%. 
Thus, we convert Ohmic loss dominated near-perfect absorber to 
near-perfect photovoltaic absorber. The power loss density 
distribution in Fig. 7(h) shows high concentration of absorbed 
power in the semiconductor layer in consistency with Fig. 7(g). 



Ohmic loss exchange at the ultraviolet frequencies. Scaling the 
absorber to the UV frequencies at Points 4 and 5 requires 
employing yet higher order surface plasmon modes with the given 
structure. However, this makes the BSDT less effective due to the 
complexity of the modes as explained below. Nevertheless, high 
absorbance of about 95% can be still achieved in the semicon- 
ductor layers mainly due to the large imaginary part of the 
permittivity of Ino.54Gao.46N. The detailed discussion for the UV 
frequencies is given in the Supplementary information. 

Discussion 

Based on simple metamaterial perfect (or near-perfect) absorber 
designs scalable to any frequency in the solar spectrum, seemingly 
contrary to common sense, we have shown that metallic layers in the 
absorbers do not necessarily constitute a problem for photovoltaic 
applications. The BSDT along with the bulk absorbance character- 
istics of the semiconductors, associated with the imaginary part of 
their complex permittivity, can allow for exchanging undesired res- 
istive losses with the useful optical absorbance in the semiconductor 
layers. Thus, Ohmic loss dominated perfect (or near-perfect) absor- 
bers can be converted into photovoltaic near-perfect absorbers with 
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Figure 7 | (a) (A/m) for a higher order plasmonic mode resonating at 672 THz (i.e., in the neighborhood of Point 3) supported by the near-perfect 
absorber with the parameters given in Table 1. Absorbance in different parts of the perfect absorber (b) operating in the neighborhood of Point 3 (legends 
not shown for clarity, same as other plots), (c) when the strip thickness in the absorber increases to fj = 40 nm, and (d) when both the strip and the ground 
plate have = tg= 60 nm thickness. All the other parameters are kept fixed and the same as in Table 1. Absorbance in different parts of the absorber 
achieved by the BSDT process described in (b)-(d) followed by the gradual incorporation of the true complex permittivity such that (e) Q(£) = 0.1, (f) 
3(e) =0.5, and (g) Q{e) =Sj(£3) = 1.61. The total absorbance in plain Iuq 54Gao.46N layer of the same thickness (i.e., d = 31 nm) is also shown for 
comparison, (h) Power loss density (W/m') distribution for the near-perfect absorber in (g) for the third selected frequency of 695 THz in the 
Ino.54Gao.46N data (i.e.. Point 3 in Table 1). Significant portion of the power (i.e., over 85%) is absorbed inside the Ino.54Gao.46N layer. 



the advantage of harvesting full potential of light management cap- 
ability offered by perfect absorber designs. It is also worth mention- 
ing that our results confirm the recent work in Ref. 34. 

The BSDT provides a simple guiding mechanism to control the 
impedance prudently in a desired fashion. If there is no such a guid- 
ing mechanism such as BSDT, with a capability alone to simply tune 
the impedance, it is difficult to control and predict in an effective 
manner where the absorbance is localized in the structure. In fact, 
this is the same reason why the BSDT as applied to some of the higher 
order modes (see, for example. Fig. S9, Supplementary information) 
is less effective; because the impedance is only partially controlled by 
the BSDT due to more complicated mode profiles of the higher order 
modes which are not fully accounted by the BSDT. For simplicity, the 
way we have applied the BSDT considers only the surface waves on 
the metal-spacer interface. However, the BSDT can be applied more 
effectively by considering all the existing surface waves constituting 
the mode. 

The perfect absorption in earlier reports^^ on metamaterial absor- 
bers was qualitatively attributed to matching the impedance 

z = <^ (j^'-eff j ^iff^ °f bulk metamaterial to that of vacuum where 

fiejf and e^jf are the bulk effective magnetic permeability and bulk 
effective electric permittivity of the metamaterial, respectively. 



However, this viewpoint has not received wide acceptance for the 
conceptual reasons as discussed below. Instead, the metamaterial 
perfect absorbers have been described more meaningfully by several 
different models such as single-resonance'^'"*^, perfectly impedance- 
matched sheet''', and equivalent circuit'^'''"''' models. In these mod- 
els, the metamaterial absorbers are defined directly by their effective 
impedance to avoid the ambiguity in defining /i^j- and since the 
metamaterial perfect absorber cannot be strictly considered homo- 
geneous bulk media. To resolve the ambiguity, in a different modeP" 
the whole metamaterial was considered as a three-layer weakly 
coupled system consisting of two metamaterials separated by a 
dielectric spacer and an interference theory''"'" was used to explain 
the perfect absorption. Non-zero order Bragg waves supported by the 
metamaterial perfect absorbers were treated in a grating theory" 
where the periodicity of the metamaterial was rigorously taken into 
account. 

The metamaterial perfect absorber structures with geometry sim- 
ilar to ours have been modeled under a single-resonance approxi- 
mation in Refs. 15 and 39. In this model, an effective surface 
impedance of the metamaterial structure is calculated. The perfect 
absorption occurs when the effective surface impedance of the meta- 
material structure matches to that of vacuum. Single-resonance 
model can be readily applied to our metamaterial absorbers 
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operating in the frequency spectrum between Points 1 and 2 since we 
use exactly the same magnetic dipole mode as in Refs. 15 and 39. 
Therefore, the single-resonance model was used in the calculation of 
the effective surface impedance shown in Fig. 4(b) to explain the 
effect of BSDT on localization of absorption in different parts of 
the metamaterial absorber. Although it is not necessary for the scope 
of the paper, the metamaterial absorbers operating around Points 3- 
5 can be also described by an appropriate model considering the 
underlying resonant modes. 

BSDT, as applied here, becomes especially effective in the large 
part of the solar spectrum below UV frequencies. Around UV fre- 
quencies, bulk absorbance characteristics of the semiconductor dom- 
inate the change in absorption localization possible with BSDT alone. 
Based on experimentally measured Ino.54Gao.46N data and a "con- 
servative" Drude model for gold^', we have shown that between 
75%-95% optical absorbance can be achieved in the semiconductor 
layers of the converted absorbers depending on the targeted fre- 
quency in the solar spectrum. However, considering more realistic 
permittivity data for gold including the size effects may result in 
slightly less shift of the absorption localization to the semiconductor 
layers than predicted here especially at small solar wavelengths. 

It is desired that optical enhancement in absorbance in the active 
region of the solar cell is maximized. This can be best achieved by 
increasing absorbance in the entire solar spectrum. Thus, the larger 
the operating bandwidth of the metamaterial photovoltaic absorber, 
the more efficient will be the solar cell. However, even a narrowband 
absorber operating around a specific frequency in the solar spectrum 
can provide overall optical enhancement in the solar spectrum com- 
pared to reference solar cell (see, for example, Ref 27); because the 
optical enhancement in absorbance achieved near the operating fre- 
quency (or frequencies if multiple resonances of a single absorber is 
considered) of the absorber is usually sufficient to obtain important 
overall enhancement in the entire solar spectrum. This is true even if 
the metamaterial structure does not perform as a perfect absorber. 
Here, we purposely tune the structure as a near-perfect photovoltaic 
absorber using BSDT. Therefore, overall optical enhancement 
should be more optimal for the near-perfect photovoltaic metama- 
terial absorbers guided by the BSDT. 

We should also note that solar cell architectures where the metals 
are completely buried inside the semiconductor can provide more 
efficient conversion process, because this increases the interaction 
surface between the metal and the semiconductor and enable more 
surface fields to interact with the electrons in the semiconductor. 
This is likely for In^Gai.^N PV consisting of multiple junctions, 
which would be necessary to obtain the highest potential efficiencies 
for this material system. This becomes especially important when 
more complex higher order surface plasmons are exploited. 
Although we chose gold and In^Gai.xN material system as an 
example to illustrate the concept of absorption localization change 
here, we can use any material systems including metals (silver, cop- 
per, aluminum, etc.), heavily doped semiconductors, nitrides'", oxi- 
jjgj44,45^ amorphous silicon, graphene and other organic materials 
provided that they have compatible fabrication process and engi- 
neered to support surface plasmons. Our approach can be also 
extended to broadband absorbers""", super absorbers'', photodec- 
tors, metamaterials^*'^^'''^, plasmonic and other metal and semi- 
conductor based optical devices where resistive losses and power 
consumption are important pertaining to the device performance. 
In addition, this approach can be applied to optimize the optical 
enhancement of solar PV devices of any arbitrary incoming spec- 
trum. For example, recent work has shown that the spectral albedo 
could play a significant role in the optimization of PV systems 
designed in snowy regions'". Following this work metamaterials 
can be designed to tailor PV devices for specific environmental con- 
ditions resulting in non-standard solar spectra (e.g. heavy pollution 
of a high UV index). 



Methods 

Simulation. All the structures shown here are designed and simulated using fully- 
vectorial finite element based commercial software package COMSOL. Gold is 
modeled by conservative Drude model using plasma frequency — 2100 THz, and 
collision frequency/^ — 19 THz^^. Due to the periodic arrangement of the plasmonic 
nanostructures, periodic boundary conditions are used in the simulations for the 
vertical boundaries and ports are used for the top and bottom boundaries. The ports 
were placed in air and sufficiently away from the metamaterial absorbers to avoid the 
evanescent fields. The distance between the input port {output port) and the top 
(bottom) surface of the metal strip (ground plate) was chosen 140 nm for the Points 
1-3 and 200 nm for the Points 4 and 5. Normally incident plane wave illumination 
with horizontal polarization (i.e., parallel to the metal-spacer interface) and 1 W 
input power, constant over the frequency spectrum, was used for the excitation. The 
power absorbed in specific regions of the metamaterial absorber was calculated by 

integrating the time-averaged power loss density [i.e., - 3?(E.J*) where E and J are 

complex electric field intensity and current density, respectively] over the specified 
regions. We also verified the results independently from power flow and s- 
parameters. The absorbance was calculated by normalizing the absorbed power in 
specific parts of the absorber with respect to input constant power. We did not use any 
reference solar spectral irradiance, since no specific solar cell architecture is 
considered here. Near Points 4 and 5, the diffraction has contribution to the 
calculated absorbance spectra. However, this does not invalidate the detailed 
qualitative discussion in the Supplementary information on absorbance localization 
but provides useful insight about the role of BSDT on complex high order modes and 
how these modes distribute the power inside the structure. 
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